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Although Godel's incompleteness theorem made mathematician recognize 
that no axiomatic system could completely prove its correctness and that 
there is an eternal hole between our knowledge and the world, and in spite 
of the work of Poincare of about 100 years ago and the further development 
of the theory of chaos, the dream of man to conquer nature and to know 
everything about nature refuse to die away. Physicists continue this ambition 
in working so far on the approaches based on the hypothesis to completely or 
approximately know the systems of interest. In this paper, however, I review 
the recent development of a different approach, a statistical theory based upon 
the notion of incomplete information. Incomplete information means that, 
with complex systems whose interactions cannot be completely written in its 
hamiltonian or whose equation of motion does not have exact solution, the 
information needed to specify the systems is not completely accessible to us. 
This consideration leads to generalized statistical mechanics characterized by 
an incompleteness parameter u which equals unity when information is com- 
plete. The mathematical and physical bases of the information incompleteness 
are discussed. 

The application of the concomitant incomplete fermion statistics to cor- 
related electron systems is reviewed. By comparison with some numerical 
results for correlated electron systems, it is concluded that, among several 
other generalizations of Fermi-Dirac distribution, only the incomplete one is 
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suitable for describing this kind of systems. The extensive incomplete fermion 
distribution n = l/{exp[u(e — ef)/ksT] + 1} gives very good description of 
weakly correlated electrons with about 0.003 < oj < 1, the normalization in- 
dex in J^iPf = 1 where pi is probability distribution. On the other hand, the 
nonextensive fermion distribution, n = 1/{[1+(uj— l)(e— ef)/kBT] u ^ u ~^ + l}, 
does not show weak correlation behaviors of electrons and is only suitable to 
describe strong correlated heavy fermion systems showing strong increase of 
Fermi momentum with increasing correlations for < u < 1. 



I. INTRODUCTION 

As the study of complexity advanced, scientists have realized that chaotic and fractal 
behaviors were ubiquitous in nature and the simple phenomena described by deterministic or 
quasi-deterministic JIJ physical sciences considering only simple interactions or predictable 
linear behaviors were only a few special or accidental cases. It was also realized that patching 
up was fundamentally useless within the conventional physics theories that break down once 
applied to complex systems having long range interactions or showing nonlinear behavior 
related to chaotic or fractal phase space structure. Generalization of these theories would 
be necessary. Driven by the increasing knowledge about chaos and fractals, the attempt of 
generalization has been rapidly focused on the problems relative to information and statistics 
theory [@-§|. The development of the nonextensive statistical mechanics (NSM) ||7|J^,p!0|-[T2 
among others | 13f , is a good example of this tendency in physics. 



Though considered by some to have a weak point due to the lack of clear physical significa- 
tions of its generalization parameter q, the probability distributions of NSM has been proved 
to be surprisingly useful for describing complex systems having long term interactions or 
correlations for which Boltzmann-Gibbs statistics (BGS) is no more valid. NSM generalizes 
BGS with a distribution function called g-exponential given by exp q (x) = [1 + (1 — q)x\ 1 ^ l ^ q ' > . 
The latter is the inverse function of a generalized logarithm hi q (x) = xl 1 l~ 1 which can be 
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used as a generalization of Hartley logarithmic information measure to obtain the g-entropy 
S q = ~k t^ q Pl (l G R ) BUI Proposed by Tsallis [0. When q = 1, These above two 
generalized functions become the usual ones and the g-entropy becomes Shannon one.[] 

In the present paper, I will review our recent efforts to find consistent foundation for NSM 
distribution functions and to give satisfactory answers to some fundamental questions. These 
efforts are based on a notion which is both new and old : incomplete information §HJ. New 
because scientists always claimed, in constructing physics theories, that their theories contain 
all necessary information for specifying the systems under consideration. This is the case of 
all the conventional physical theories : from Newtonian to quantum physics, in passing by 
Einstein, Boltzmann and Shannon (certainly, a theory containing only partial information 
about the system of interest is a little bit discouraging). Old because since the discovery 
of, e.g., irrational numbers, mathematicians know that, within arithmetical system, they 
loss some information about the world and that one could not know everything with infinite 
precision. In 1931, Godel shown [f§-|6|] that mathematics system (or any axiomatic system) 
is incomplete in the sense that within any such axiomatic system there is never sufficient 
information to prove all possible statements of the theory . If a non negligible amount of 
information is not accessible to us, BGS theory has to be modified. Incomplete information 
theory is a kind of modification (generalization) of BGS suggested by this consideration as 
well as by some difficulties encountered within NSM in the last decade [El 151. 



From now on, the parameter q will be replaced by uj and Tsallis entropy by S u = —k — 1 ■ 

The above generalized functions will be called ^-exponential and w-logarithm. I make this replace- 
ment for the simple reason that, though it often gives similar forms of functions as q, uj defined in 
the framework of the theory I review here does not have the same physical content as the parameter 
q in Tsallis version of NSM. So I prefer to use u to avoid confusions. By definition, u> has clear 
physical meaning as the reader will find in this paper. 



3 



II. COMPLETE INFORMATION ASSUMPTION 



In this section, I will briefly review the well known information theory founded by Shan- 



non et al [ 16 1 • It should be remember that information about a real system is not our 
knowledge about it. It is our ignorance. The ignorance of something to which we may have 
access. A mail address, as a state of physical system, may be an information if we do not 
know it. More we know about a system, less there is information in its description. So in 
a deterministic theory (e.g., classical mechanics), information is null. In statistical theory, 
there is information because we ignore something so that we are not sure of the exact state 
at any given moment of the system under consideration. So information can be related to 
the uncertainty due to the ignorance or to the probability of finding the system at different 
states. It should be noted that, as mentioned above, up to now, we always suppose that the 
information we address in any statistical theory is complete or completely accessible. That 
is if we obtain it, we can answer all questions which can be asked about the system. This 
certainty is reflected by the following postulate : 

i=i 

where v must be the number of all the possible states of the system under consideration. As 
a result, the arithmetic average of £ is given by x = Y%=iPi x h 

By some analysis of the information properties, it is supposed PJT6|] that the information 
is given by the well know Hartley formula ln(iV) |l7j needed to specify N elements, or by 
ln(l/pj), the information needed to specify that an element will be found at the state i. If 
we perfectly know all the v possible states, then the complete information measure I is given 
by averaging all ln(l/pj) : 

I = 5>ln(l/ Pi ). (2) 

i=i 

It should be emphasize that the above definition of information or entropy needs the harsh 
condition that the interactions in the system of interest are of short range or limited between 
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the walls of the containers of subsystems which are consequently independent of each other. 
To see this, it suffices to consider the assumption of information additivity, i.e., for a system 
C containing two subsystems A and B, it is supposed 1(C) = 1(A) + 1(B). This additivity 
is valid if and only if the information 1(C) needed in order to specify simultaneously A and 
B is given by hi[N(A)N(B)] where N(A) and N(B) are respectively the number of elements 
in A and B. This is as if we had a system C containing N(A)N(B) elements. This result 
needs that the states of the elements of A do not depend on the states of B. In other words, 
these is no interactions between the elements of A and those of B. There may be interactions 
between the elements on the walls of the containers of A and B, but most of the elements 
inside A and B must be independent. This is a case of short range interaction where we 
have not only additive information or entropy, but also additive energy and other extensive 
thermodynamic variables. 

I would like to recall in passing here that the total information \n[N(A)N(B)] implies 

p ij (C)= Pi (A)p j (B) (3) 

where Pij(C) is the probability that the composite system C is at the product state ij when 
A is at the state i with probability Pi(A) and B at j with Pj(B). Eq.(|3|) symbolizes the 
independence of the noninteracting subsystems having additive physical quantities. But for 
interacting subsystems, it symbolizes totally different physical reality. This product law has 
been widely employed and discussed in the last decade in connection with equilibrium and 
many body problems |T8|-p4| and caused much confusion within NSM because it paradoxi- 
cally independence of subsystems and additive energy for nonextensive interacting systems. 
Very recently, we shown that Eq.(|]) was nothing but the consequence of the existence of 
thermodynamic equilibrium in interacting systems described by u -entropy and did not need 
independence of the subsystems. This conclusion allows to exactly define equilibrium param- 
eters such as temperature, pressure and chemical potential for nonextensive systems and to 
obtain the exact one body quantum distributions [fZT|-pi|]. 

According to above discussions, we can say that, if there are long range interactions 
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between A and B, the information about C will be different from \n[N(A)N(B)] because the 
elements are correlated and can no more occupy their states independently. According to the 
nature of the correlation, there may be more or less information than in the noninteracting 
case. In general, we should write 1(C) = 1(A) + 1(B) + f [1(A), 1(B)], a case treated by 
NSM. Now Eq.(|3]) becomes questionable, yet it is a crucial relationship for any statistical 
mechanics, for it's applications to many-body systems and it's thermodynamics connection. 
The reader will find detailed discussions on this issue below. 



III. COMPLEXITY AND MATHEMATICS 

Certainly, complete information is possible whenever all possible states are well known 
so that we can count them to carry out the calculation of probability and information. In 
physics, this requires that we can find the exact hamiltonian and also the exact solutions 
of the equation of motion to know all the possible states and to obtain the exact values 
of physical quantities dependent on the hamiltonian. The reader will see that these two 
"exact" conditions of complete information are almost impossible to satisfy. 

Let us begin by asking some questions about the mathematical basis of physical theory. 

What is the A basic field of mathematics is the classical arithmetic. From the episte- 
mological point of view, arithmetic is a theory based on a model of world resulted from the 
direct intuition of human beings. This is a simple model for fragmented world containing 
only isolated, distinct and independent parts. So you have 1 = 1, 1 + 1 = 2 and a series 
of rules, theorems and generalizations. No matter how complicated are the immense math- 
ematical constructions developed from arithmetic, their validity is always limited by these 
initial conditions imposed by the crude data of our senses and direct intuition. Indeed, our 
senses, luckily, have the capability of filtering the complex world into separated and dis- 
cernible parts. If not, scientific knowledge would be impossible. But these harsh constraints 
imposed by this filtration, as claimed by Poincare ||25| , should not be forgotten. We have 



to ask the following question : how far he can go with the concepts formed through the 
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filtration in the real messy world or complex systems including interacting, entangled and 
overlapped parts, especially when the interactions can no more be neglected. 

So in some sense, it can be said that mathematics is an approximate theory contain- 
ing finite amount of information about the world which is surely incomplete because some 
information is lost by our senses through the formation of the axioms. Any formation of 
axiomatic systems is necessarily made through a kind of filtration of the world. The results 
of the filtration are not wrong, but they are only partially true. Something about the con- 
nection of different parts of the world is rejected by the filtration. In my opinion, this is 
why axiomatic systems, as stated by the incompleteness theorem of Godel, inevitably fail to 
prove some statements, especially those about their axioms. There is no enough information 
for that. The missing information is just what rejected by the formation of axioms. 

A mathematician is rather interested by the coherence of his logical systems based on 
axioms. He may put aside the missing information and work within the logical systems 
without being connected to physical reality. But for a physicist, the connection of his theory 
to the outside world is the most important thing he mind. He possibly ask : My physical 
theory is in fact an application of a incomplete mathematical theory. If the information I 
am handling is not complete, how can I apply it to the world whose description probably 
needs more information? 

In what follows, we will try to answer this question in recognizing that the incompleteness 
of all axiomatic systems discovered by Godel has put an end to the ambition of establishing 
physical theories containing or capable of treating complete information about any system 
in the world. In this sense, any physics theory is incomplete by definition. This is the 
very reason for the introduction of "incomplete information" into statistical physics. This 
introduction needs in addition other considerations I am presenting below. 
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IV. COMPLEXITY AND INCOMPLETE INFORMATION 



Now let us look at the information problem from the physical viewpoint. I will try to 
show that, due to the omnipresent complexity in the world, we cannot have access to all the 
necessary information for complete description of a system. Here "complexity" means that 
the systems show nonlinear behaviors which are extremely sensible to initial conditions and 
unpredictable. This is the famous chaos observed almost everywhere in the world 

A complex system is not necessarily a complicated system with a large number of free- 
doms. A one dimensional oscillator with well known nonlinear interaction (with potential 
oc x 4 , for example) or a three body system with gravitation (oc 1/r) can behave chaotically. 
These two cases are just very good examples of the impossibility of the two "exact" con- 
ditions of complete information mentioned above. In the case of the three body problem, 
we know (at least we believe that we know) the exact interaction of the system (Newtonian 
gravitation). But Poincare showed that the exact and predictable solution of the equation 
of motion was not possible . There are in fact infinite number of periodic and aperiodic 
solutions. The movement is chaotic and unpredictable and the attractors of the chaotic 
structures formed by the trajectories in phase space are fractal. This means that we never 
know all possible states of the system and that complete information treatment becomes 
impossible. We even have to redefine probability distribution in order to calculate it in 
chaotic or fractal phase space. 

Above conclusion is for hamiltonian systems whose interactions is a priori well known. 
When the hamiltonian cannot be exactly written, the situation is more complicated. Even 
the exact and predictable solutions of equation of motion are not complete due to the 
incomplete hamiltonian. This may happen if, for a isolated closed system, the interactions 
are too complex to be written, or, for a system with simple interactions, the effects of 
the external perturbations are not negligible. Sometimes negligible perturbations may have 
drastic consequences if the system is sensitive to initial conditions. In this sense, the omission 
of small interactions may make enormous information unaccessible to the theory. This 
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incompleteness due to neglected interactions simply adds to the incompleteness mentioned 
previously. 

In any case, complete information description of complex systems is only a science fiction. 
Although we cannot say that all these systems have chaotic or fractal nature, a common 
feature of them is that a part of their phase space is unknown so that complete and exhaustive 
exploit of the phase space is impossible. The calculable information is inevitably limited 
by this incompleteness of knowledge. That is evident. The treatments of these systems 
based on the assumption of complete information and probability distribution are not well 
founded. They are legitimate only when unaccessible part of the information is negligible 
with respect to the accessible information and to the desired precision of observation or 
theoretical description. 

In what follows, we will try to introduce the notion of incompleteness of information into 
physics through statistical method. It was with this method that man began to overcome 
the obstacle of his limited knowledge in supposing, on the basis of Newtonian or quantum 
mechanics, that the missing knowledge (information) is mathematically accessible or, equiv- 
alently, that the calculated probability must sum to one. Now if we say that we cannot have 
access to every information we need or to every point of the phase space, a serious impact 
on the normalization of probability, the very first stone in the construction of statistics, will 
be inevitable. 

V. CHAOS AND INCOMPLETE PROBABILITY DISTRIBUTIONS 

A. Incomplete normalization 

What can we do for probability and information calculation if we do not know how 
many states the system of interest has? When we deal with a chaotic system having fractal 
attractor in phase space ||, it is as if we toss a coin which often comes down, neither tails 
nor heads, but standing on the side without, in addition, being observed. All calculations 
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based on Eq.(|l]) with v = 2 would lead to aberrant results because we have now Yn=iPi — 
Q / 1. In this case, Pi is referred to as incomplete distribution || and Q is a constant 
depending on and characterizing the incompleteness of the system and provides a possible 
key to introduce incompleteness of information into physics theory. It should be supposed 
Q = 1 if information is complete. 

The philosophy of incomplete information theory we developed is to keep the methods of 
classical complete probability theory for incomplete information or probability distribution 
by introducing empirical parameters in order to characterize the incompleteness. This is just 
the same methodology as in the theory of chaos or fractals introducing fractal dimension 
to characterize the structures of space time. In this sense, we can refer to the parameter u 
introduced below as incompleteness parameter. 

First of all, we need a "normalization" for incomplete distribution p, in order to take 
advantage of the conventional probability theory. This is an occasion to introduce a 
parametrization function and to write 

E^fe) = 1 (4) 

i 

which can be called generalized or incomplete normalization. F u should depend on the nature 
of the system and become identity function whenever information is supposed complete 
(Q = 1). The arithmetic average should now be given by x = J2i F^ip^Xi. F^ can be 
determined if the information measure and the distribution law are given. For example, 
with Hartley information measure and exponential distribution, F^ can be showed to be 
identity function |J. In general, by entropy maximization through the functional 

8[Y,FMi{Pi) + PY. F M x i\ = ( 5 ) 

i i 

we get : 

d In F w fa) dl/dpi , 

dpi ~I + (3f- l {pi) U 

or F u (pi) = C exp[/ rr^T% \ dpi\ where {3 is the multiplier of Lagrange connected to expec- 

l^~PJuj \Pi) 

tation, I(pi) is the information measure, Pi = f^Xi) the distribution function depending on 
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the parameter oj, C the normalization constant of F u . 

B. Incomplete normalization of NSM 

In my previous papers in order to find coherent foundation for cu-exponential dis- 

tribution on the basis of cu-logarithm information measure, F^ipi) = pf was postulated. So 
that 

E# = 1 - (7) 

i 

In what follows, I will try to show that the conjecture of power law incomplete normalization 
in the previous section is inevitable in a chaotic or fractal space time. 

For the sake of simplicity, let us consider a phase space in which the trajectory of a 
chaotic system forms a simple self-similar fractal structure, say, Sierpinski carpet (Figure 
1). This means that the state point of the system can be found only on the black rectangular 
segments whose number is Wk = 8 k at k th iteration. Hence the total surface at this stage is 
given by Sk = W^su where = lo/3 h is the surface of the segments at k th iteration and Iq the 
length of side of the square space at th iteration. If the segments do not have same surface, 
we should write Sf. = Y^i s k{i)- We suppose that the density of state is identical everywhere 
on the segments and that the distribution is microcanonical, so that the probability for the 
system to be in the i th segment may be defined as usual by pi = Sk(i)/Sk- This probability 
is obviously normalized. The problem is that, as discussed in ||, Sk is an indefinite quantity 
as k — > oo and, strictly speaking, can not be used to define exact probability definition. 
In addition, Sk is not differentiate and contains inaccessible points. Thus the probability 
defined above makes no sense. 

Alternatively, the probability may be reasonably defined on a integrable and differentiable 
support, say, the Euclidean space containing the fractal structure. To see how to do this, 
we write Sk = lo(-^k) d ~ df for identical segments or, for segments of variable size, 

E[^V' /d = 1 (8) 
i=i 
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where Sq = Iq (here d = 2 for Sierpinski carpet) a characteristic volume of the fractal 
structure embedded in a <i-dimension Euclidean space, df — is the fractal dimension, 
n = 8 the number of segments replacing a segment of the precedent iteration and m = 3 the 
scale factor of the iterations. The microcanonical probability distribution at the k th iteration 
can be defined as pi = so that Ya^ipI 1 ^ — 1 which is just Eq.(|7|) with uj = df/d. The 
conventional normalization Y^=iPi — 1 can be recovered when df = d. 

It should be noticed that, in Eq.(||), the sum over all the Wk segments at the k th iteration 
does not mean the sum over all possible states of the system under consideration. This is 
because that the segment surface Sfc(i) does not represent the real number of state points on 
the segment which, as expected for any self-similar structure, evolves with k just as S^. So 
at any given order k, the complete summation over all possible segments is not a complete 
summation over all possible states. But in any case, whatever is k, Eq.(||) and Y^iPt — 1 
always holds for u = df/d. 

In this simple case with self-similar fractal structure, the incompleteness of the normal- 
ization Eq.(|7]) is measured by the parameter u = df/d. If df > d, there are more state 
points than Wk, the number of accessible states at given k. If df < d, the number of ac- 
cessible states is less than Wk- When df = d, the summation is complete at any order k, 
corresponding to complete information calculation. 

VI. INCOMPLETENESS PARAMETER u 

Here I will discuss in a detailed way the incompleteness parameter u and its physical 
meanings. Incomplete statistics gives to the empirical parameter u a clear physical signifi- 
cation : measure of the incompleteness of information or of chaos. Let us illustrate this by 
the simple case of self-similar fractal phase space with segments of equal size. 
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A. u and phase space expansion 

As discussed in the case of chaotic phase space, u = In n/d In m gives a measure of the 
incompleteness of the state counting in the <i-dimension phase space, u = 1 means df — d 
or n = m d . In other word, at the k th iteration, a segment of volume Sk is completely covered 
(replaced) by n segments of volume Sk+i = Sk/m d . So the summation over all segments 
is equivalent to the sum over all possible states, making it possible to calculate complete 
information. 

When u > 1 (or u < 1), n > m d (or n < m d ) and Sk is replaced by n segments whose 
total volume is more (or less) than s&. So there is expansion (or negative expansion) of state 
volume when we refine the phase space scale. An estimation of this expansion at each scale 
refinement can be given by the ratio r = nSk + 1 ~ Sk = ™ _ 1 = — 1 = (hj — 1) I^-J — — — . 

r describes how much unaccessible states increase at each step of the iteration or of the 
refinement of phase space. The physical content of u is clear if we note that uo > 1 and 
uj < 1 correspond to an expansion (r > 0) and a negative expansion (r < 0), respectively, of 
the the state volume at each step of the iteration. When u = 0, we have df — and n — 1, 
leading to r = — 1. The iterate condition n > 1 means lo > 0, as proposed in references 
||. uj < is impossible since it means df<0orn<l which obviously makes no sense. We 
can also write : u — 1 = ln(r + l)/ln(m d ) = ln(nsfc + i/sfc)/ln(m d ), which implies that it is 
the difference u — 1 which is a direct measure of the state space expansion through the scale 
refinement. 

B. u and information growth 

The expansion of the state volume of a system in its phase space during the scale refine- 
ment should be interpreted as follows : the extra state points A = nsk+i — Sk acquired at 
(k+ l) th order iterate with respect to k th order are just the number of unaccessible states at 
k th order. A > (or A < 0) means that we have counted less (or more) states at k th order 
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than we should have done. A contains the accessible information gain (AIG) through the 
(k + l) th iterate. 

To illustrate the relation between this "hidden information" and the parameter u, let us 
first consider the Hartley logarithm information in the simple case where the distribution is 
microcanonical and scale-invariant |f2B[. At the iterate of order k, the average information 
contained on Sk is given by Ik = J Sk p u \a.{l/p)ds. At k + 1 order, Ik+\ = J nSk+1 P^ hi(l/p)ds. 
Hence AIG is just AJ = I k+ i — h = I(ns k+1 -s k ) P^ hi(l/p)<is = O/A, where 07 = p w ln(l/p) is 
the information density or the average information carried by each state. The relative AIG 
is given by AJ/J^ = r = (1 — uj) ^^^- ~ 1 which is independent of scale but dependent on 
scale changes. For given scaling factor m, the magnitude of AI or r increases with increasing 
difference |1 — oj\. The sign of r (or AIG) was discussed earlier. For given uj, \AI\ increases 
with decreasing scaling. For uj = 1 or m = 1, there is no information gain, corresponding to 
the case of complete information. 

According to the relationship u = df/d and the above discussions, it can be concluded 
that the incompleteness parameter u may be considered as a measure of chaos. Certainly 
this is a conclusion on the basis of simple models and the relation between uj and the degree 
of chaos or fractal may be more complicated with more complex chaos and fractals, but it 
is consequent to say that more a system is chaotic, more its information is incomplete and 
more uj is different from unity. 



VII. NON ADDITIVE INCOMPLETE DISTRIBUTIONS 

To get the nonextensive distribution in tu-exponential as mentioned above, we can max- 
imize the entropy S w = — k ^- JiPl 1 ^ iPi |||| according to the Jaynes principle |27| with the 
constraints U = J^iPfEi and N = J^iPfNi for grand-canonical ensemble, where U is the 
internal energy, N the average particle number, E{ the energy and A, the particle number 
at the state % of the system. We obtain : 

Pi — ' \ ' 
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where Z w = ~ (1 ~ ~ ^i)]t~"- Wt = a; if x > and [x] t = otherwise. 

/3 is the inverse temperature and \x the chemical potential. This distribution function has 
been proved particularly useful for systems showing non gaussian distribution functions 
(for detailed information, see JTTJ and references there-in). Considering Eq. ([!]), the product 



probability law at thermodynamic equilibrium, the one-particle distribution from Eq.(^) can 

be rewritten as pk = — where is the energy of one particle at the state k 

and z^ = — (1 — uj)(3{ek — fj)]l~ u is the one-particle partition function. 



As shown in |[23|| , the above one-particle distribution can be recast into exponential form 
as follows 

Pk = -[1 - (1 - uj)P'e k }^[l + (1 - cu)(3'fz']T^ = L e -P'^ (10) 

Z Zj 

where j3' = l _^ 1 ~ )/3fi , n' = /i[l-(g-l)/3/i] which imply /3'// = /5yu,, v = ^^'j^f - 1 and e fc = 
^T^zjj^r £fc ^ ■ The exponential distribution Eq.(|10"D makes it possible to straightforwardly 
obtain the exact quantum distribution |23|] (EQD) given by 

where is the occupation number of the one-particle state k. "+" is for fermions and "-" 
for bosons. These distribution can be compared to the approximate quantum distributions 



(AQD) of NSM rifc = m \ , — = r given within a factorization ap- 

e k ±X [l+( 9 -l)/J(e fe - M )] — ±1 

proximation using additive energy. At first glance, EQD and AQD are not very different 
from each other if we put ui = q. But Figure 2 shows that they are two very different distri- 
butions. AQD remains approximately the same as the conventional Fermi-Dirac distribution 
for whatever q value. So its Fermi energy e/ is almost constant with changing q. On the 
contrary, EQD changes drastically with u. The Fermi energy Cf shows a strong increase 
with decreasing uj up to two times e^ of the conventional Fermi-Dirac distribution when 
u ^ 0. This ej increase has been indeed noticed through numerical calculations for strongly 
correlated heavy electrons on the basis of tight-binding Kondo lattice model |2"H| , |2D|] as shown 



in Figure 3. Increasing correlation corresponds to decreasing u from unity (zero correlation). 
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This implies that EQD based on incomplete information has its merit in the description of 
heavy electron systems. Further investigation is needed to know the connection between the 
correlation and the nonextensive parameter 1 — u. 



VIII. ADDITIVE INCOMPLETE DISTRIBUTIONS 

Although the nonextensive EQD accounts for an important aspect of correlated electrons, 

1. e., the correlation induced Fermi energy increase, another important aspect of the weak 
correlation is missing in the description of nonextensive EQD. This is the flattening of n drop 
at e/ p8|-[3l|]. That is the correlation, even at low temperature, drives electrons above e/ so 
that the n discontinuity becomes less and less sharp as the correlation increases. Curiously, 
this flattening of n discontinuity at e/ is completely absent in EQD of NSM. From Figure 

2, we see that the sharp n drop at e/ is independent of u or correlations. 

In what follows, I will present an additive incomplete statistical mechanics. It is assumed 
that the additive Hartley information measure still holds. So with respect to the conventional 
Shannon information theory and BGS, only the normalization is changed according to Eq.([7|) 
|9|Jl3||. The additive incomplete entropy is given by S = kY,T=iPi hi(l/pj). When oj — > 1, S 



is Shannon entropy, which identifies k to Boltzmann constant. 

For grand canonical ensemble, the usual entropy maximization procedure leads to 
Pi = Q-uPiEi-nNi) i z where part i t i on f unc tion is given by Z = {TX=i e-^ Ei ~^ . For 
quantum particle systems, we have 

nk = e^»-t)±r (12) 

The fermion distribution given by Eq.(|T2"|) is plotted in Figure 3 for different uo values 
in comparison with some numerical simulation results. We note that IFD reproduces well 
the numerical results for about J < 1. When coupling is stronger, a long tail in the KLM 
distributions begins to develop at high energy. At the same time, a new Fermi surface 
at k = kf + 7r/2 = 0.757T starts to appear and a sharp n drop takes place at the new 
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Fermi momentum. At J = 4, KLM distribution (x-marks) is very different from IFD (e.g. 
u = 0.0011). The solid line fitting better the J = 4 KLM distribution is given by the 
incomplete statistics version of fractional exclusion distribution (1/n— a) a (l/n — a + l) l ~ a = 
e u)P(e-e f ) p2||33| w jth 1/a — 0.85 due to the KLM occupation number smaller than 0.5 at 
low momentum k. 

IX. CONCLUSION 

Summing up, I have discussed the philosophical basis of incomplete information from 
both the viewpoints of mathematical and physics. The information we deal with in scientific 
theories can not be complete in the sense that a part of the information necessary for 
complete description of the system under consideration is not accessible to our theory or 
knowledge. This part of information is rejected from scientific knowledge by the formation of 
concepts, axioms and models. The amount of rejected information is particularly important 
for complex systems having chaotic behaviors and fractal phase space. A parameterized 
normalization Y^i Pt — 1 is proposed for this kind of systems, where u is the incompleteness 
parameter characterizing the inaccessibility of phase space points or of the information of 
the system. It also offers a measure of the degree of chaos. 

The wide drop in the fermion occupation number and the sharp cutoff of occupation 
number at e/ showing strong increase with increasing interaction can be interpreted by the 
nonextensive incomplete fermion distribution with decreasing uj value. On the other hand, it 
fails to describe weak correlation effect on electrons which is well accounted for by additive 
incomplete fermion distribution. But the additive distribution does not show the sharp 
cutoff at ej when correlation is strong. This result suggests to combine these two partially 
valid models to describe correlated electrons in a global way. Further results of this current 
work will be presented in other papers of ours. 
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Figure caption : 

FIG. 1. A simple model of fractal phase space in Sierpinski carpet (or sponge). At k f h iteration, 
the side of the squares (black or white) is = lo/2 k and their number is Wk = 8 k , Ik being 
the length of the side at at th iteration. The total surface at k f h iteration is Sj, = WkSk or 
WkSk/Sk = 1. The classical probability definition by relative frequency of visits of each point by 
the system must be modified because the total number of visits (propotional to black surface Sk 
of the carpet) is no more a finite quantity. (Construction of Sierpinski carpet. First iteration c(l) 
: removing the central square formed by the straight lines cutting each side into three segments of 
equal size. Repeat this operation on the 8 remaining squares of equal size and so on.) 

FIG. 2. Nonextensive fermion distributions of AQD and EQD of incomplete statistical mechan- 
ics. AQD distribution is only slightly different from that at q = 1 (conventional Fermi-Dirac 
distribution) even with q very different from unity. But EQD changes drastically with decreasing 
u). As lo — > 0, the occupation number tends to 1/2 for all states below ej which increases up to 2 
times ey , the conventional fermi energy at T = 0. 

FIG. 3. Comparison of additive incomplete fermion distribution (IFD, lines) with the numerical 
results (symbols) of Eder el al on the basis of Kondo lattice t — J model (KLM) for different 
coupling constant J [Phys. Rev. B, 55(1997)6109]. In my calculations, the density of electrons 
is chosen to give fe/ = 0.257T in the first Brillouin zone. We note that IFD reproduces well the 
numerical results for about J < 1. When coupling is stronger, a long tail in the KLM distributions 
begins to develop at high energy. At the same time, a new Fermi surface at k = kf + tt/2 = 0.757T 
starts to appear and a sharp n drop takes place at the new Fermi momentum. At J = 4, KLM 
distribution (x-marks) is very different from IFD (e.g. u = 0.0011). The solid line fitting better 
the J = 4 KLM distribution is given by the incomplete statistics version of fractional exclusion 
distribution (1/n - a) a (l/n - a + 1) 1 "" = e^ e ~ e s) [Yong-Shi Wu, Phys. Rev. Lett, 73(1994)922] 
with l/a = 0.85 due to the KLM occupation number smaller than 0.5 at low momentum k. 
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Figure 1 : A simple model of fractal phase space in Sierpinski carpet. At $k A th$ 
iteration, the side of the squares (black or white) is $l_k=l_0/2 A k$ and then- 
number is $W_k=8 A k$, $1 k$ being the length of the side at at $0 A {th}$ 
iteration. The total surface at $k A th$ iteration is $S_k=W_k s_k$ or $W_k 
s k/S k=l$. The classical probability definition by relative frequency of visits 
of each point by the system must be modified because the total number of visits 
(propotional to black surface $S_k$) is no more a finite quantity. 

(Construction of Sierpinski carpet. First iteration c(l) : removing the central 
square formed by the straight lines cutting each side into three segments of equal 
size. Repeat this operation on the 8 remaining squares of equal size and so on.) 
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